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Adiabatic processes driven by non-Hermitian, time-dependent Hamiltonians may be sped up 
by generalizing inverse engineering techniques based on Berry's transitionless driving algorithm 
or on dynamical invariants. We work out the basic theory and examples described by two-level 
Hamiltonians: the acceleration of rapid adiabatic passage with a decaying excited level and of the 
dynamics of a classical particle on an expanding harmonic oscillator. 
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I. INTRODUCTION 



Non-Hermitian Hamiltonians: basic formulae 
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We refer to fast time-dependent processes that repro- 
duce the effect of a slow, adiabatic driving of a quantum 
system as "shortcuts to adiabaticity" |ll-l23j . We also ap- 
ply the term to the inverse engineering methods used to 
design these processes. In the adiabatic process of refer- 
ence the external control parameters are modified slowly 
from some initial configuration to a final one. In the cor- 
responding shortcut the system is driven in a predeter- 
mined short time to a final state which reproduces in the 
instantaneous basis the initial populations, as the adia- 
batic process would do, but possibly allowing for some 
transient excitation along the way. There is nowadays 
considerable interest in these questions for fundamental 
and practical reasons. Adiabatic methods are ubiquitous 
in cold-atom and atomic-physics laboratories to manip- 
ulate and prepare atomic states in principle in a robust 
way. An obvious drawback is that the times required may 
be too long for practical applications. Moreover the ideal 
robustness may be spoiled by the accumulation of per- 
turbations and decoherence due to noise and undesired 
interactions. Studies and experiments to speed up adia- 
batic processes have been carried out for transport [I-p , 
wave splitting 0, expansions and compressions [3- 
^d\, or internal state control [2T-'23']. These studies have 
so far been performed for Hermitian Hamiltonians, but 
many systems admit an effective non-Hermitian descrip- 
tion. In this paper we put forward shortcuts to adiabatic- 
ity techniques for non-Hermitian Hamiltonians. Specifi- 
cally we shall generalize the inverse engineering method 
proposed by Berry and the one based on dynamical 
invariants While these methods are intimately con- 
nected as shown in [23*1 and may in fact be considered 
potentially equivalent, in standard applications they are 
used in different ways and provide different answers so 
we shall consider them separately here. As study cases 
we shall discuss a two-level decaying atom and the mo- 
tion of a classical particle in a harmonic oscillator with 
time-dependent frequency. 



Non-Hermitian Hamiltonians typically describe sub- 
systems of a larger system 0|. We shall first review 
a basic set of relations and notation [23|. We shall 
assume a non-Hermitian time-dependent Hamiltonian 
Ho^t) with N non-degenerate right eigenstates {\n{t))}, 



n = 1,2..., A^, 



Ho{t)\n{t)) = En{t)\n{t)), 
and biorthogonal partners {\h{t))}, 

Hiit)\h{t)) ^ E:{t)\m), 



(1) 



(2) 



where the star means "complex conjugate" and the dag- 
ger denotes the adjoint operator. They satisfy 

{h{t)\m{t)) ^ Snm (3) 

and the closure relations 

J2\nit)){nit)\=J2\nit)){hit)\^l. (4) 



(5) 
(6) 



{fi{t)\ is the left eigenvector of Ho{t), 

{n{t)\Ho{t) = {n{t}\En{t), 
and {n{t)\ the left eigenvector of i?o(i), 

{nit)\Hl,it) = {nit)\E:it). 
We can thus write the Hamiltonian and its adjoint as 
Ho{t) = ^|n(t))i?„(t)(n(i)|, 

n 

Hlit) = J2\^(^))^nmn{t)\. (7) 

n 

The time-dependent Schrodinger equations for a generic 
state and for its biorthogonal partner satis- 

fying = 1 are 



(8) 
(9) 
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II. TRANSITIONLESS DRIVING ALGORITHM 



as 



In [21| M. V. Berry proposed a method to design a 
Hermitian Hamiltonian H{t) for which the approximate 
adiabatic dynamics driven by the Hermitian Hamihonian 
HQ{t) becomes exact. We shah generahze this method 
for non-Hermitian Hamihonians. First we need the 
adiabatic approximation when -ffo(0 is non-Hermitian 
[ISIlll- A general time-dependent state is a hnear 

combination of instantaneous eigenvectors \n(t)) of i?o(0 
with time-dependent coefficients. Similarly \'^{t)) is a lin- 
ear combination of instantaneous eigenvectors \n{t)) of 
HQ{t). In the adiabatic approximation we assume that 
only one of these eigenvectors is populated. To determine 
the corresponding phase factor we insert 



into Eqs. (O and ©. Thus we have 

^P^\n{t)) + \^Mt)) = ^^\nit)), 



ifin\n{t)) + \dtn{t)) 



ih 



h{t)), 



(10) 
(11) 



(12) 
(13) 



where the dot denotes the derivative with respect to time. 
Multiplying Eq. ^ by {h{t)\ and Eq. ^ by {n{t)\, 
taking into account Eqs. ((T|) and ([2]), and integrating, we 
find 



Pnit) 



h 

-Kin 



+ i{n{t')\dt'n{t')) 
+ i{n{t')\dt'h{t')) 



dt', (14) 
dt', (15) 



where the initial phases are set to zero. As {n{t)\dtn)* = 
{dtn\n{t)) and, from Eq. ©, {dth\n{t)) = -{h{t)\dtn), 
we have that /3„ = /3*. 

As in [2l|, we now impose that all |(p„(i)) satisfy ex- 
actly the Schrodinger equation for a yet unknown H{t), 



Similarly, 



(16) 



(17) 



The states \ipn(t)) and |(^„(t)) can be written in terms of 
the corresponding evolution operators U{t) and U{t), 



\Vn{t)) 



U{t)\n{0)), 
U{t)\h{0))- 



The Hamiltonian H{t) can be found from 
ihdtU{t) = H{t)U{t), 



(18) 



(19) 



H{t) = ihdtU{t)U\t), (20) 

since [131 

UHt)U{t) = lop. (21) 
The evolution operators can be written as 
U{t) = ^e^^"W|n(t))(n(0)|, 

n 

Uit) = ^e^^"(*)|n(i))(n(0)l. (22) 



(23) 



Using now Eq. (|20)) . 
where 

iJi(i) = ihY,[\dtn{t)){h{t)\ 

n 

- {n{t)\dtnmn{t)){nm- (24) 

H{t) drives the system along the adiabatic paths defined 
by/fo(i). 

As noted in [21] and [23] this Hamiltonian is not 
unique. For a given set \n{t)) the same final populations 
are found by choosing different phases. Let us rewrite 
\'^n{t)) and \Lpn{t)) in terms of arbitrary phases, ^n(i) 
and £,n{t), which we now consider manipulable functions 
obeying ^„(i) = il{t) so that {(pn{t)\fn{t)) = 1, 

Wn{t)) = e*«"(*V(t)), |^„(i)) = e'«"(*)|n(t)). (25) 

We assume fn(0) = Cn(0) = and define the new evolu- 
tion operators 

m) - E^^^"^*V(0)("(o)i, 

m) = (26) 

n 

From Eq. f|2Qp . the corresponding Hamiltonian becomes 
HS) = "^E H-t))Ut){n{t)\ + |9tn(i))(n(i)l. 

n n 

(27) 



III. TRANSITIONLESS DRIVING ALGORITHM 
APPLIED TO A DECAYING TWO-LEVEL ATOM 

A. Hi{t) applied to a decaying two-level atom 

As an example of the approach of the previous sec- 
tion we shall speed up adiabatic processes in a two-level 
atom with spontaneous decay. If the decayed atom es- 
capes from the trap by recoil, a Hamiltonian (rather than 
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master equation) description is enough. We shall also as- 
sume a semiclassical treatment of the interaction between 
a laser electric field linearly polarized in cc-direction, and 
a decay rate (inverse life-time) F from the excited state. 

Applying the electric dipole approximation, a laser- 
adapted interaction picture and the rotating wave ap- 
proximation, the Hamiltonian is 



(28) 



in the atomic basis |1) 



The 



detuning from the atomic transition frequency ujq is 
A{t) = ujQ ~ i-Oi{t), where uJi{t) is the instantaneous field 
frequency. We assume a slowly varying pulse envelope 
so that the Rabi frequency riij(i), assumed real, depends 
on time. In the example below we shall take F as a con- 
stant although, in a general case, it could also depend on 
time, F — T{t), as an effective decay rate controlled by 
further interactions, see e.g. The eigenvalues of this 
Hamiltonian are 

E±it) = ^ |-zF ± y/-[F + 2zA(t)]2 + mlit)^ , (29) 

and the normalized eigenstates are 

|x+(t))=sin(|)|l)+cos(|)|2), 
|X-(t))=cos(|)|l)-sin(|)|2), (30) 

where the mixing angle a = a(t) is complex and defined 
as 



tana 



R 



A-iT 



The adjoint of Hao{t) is 



h ( -Ait) nn{t) 
2 I flRit) A{t) + iT 



(31) 



(32) 



with eigenvalues E^. (t) and normalized eigenstates 

|x+(t))=sin(^^)|l)+cos(^^) |2), 

IX-W)-cos(^^)|l)-sin(^^)|2). (33) 

Note that the coefficients are complex conjugate of those 
in Eq. ([50]) because Hao{t) is equal to its transpose [23 |. 
For this system Eq. ([M)) takes the form 

Haiit) = zh[\dtx+{t)){x+{t)\ 

- {x+{mx+mx+{t)){x+{t)\ 

+ \dtx-it)){x-{t)\ 

- {xAmx-mx-{t)){x-im, m 



where, according to Eqs. pOI) and 

{x±(t)\dtx±it)) = 0, 



{xT(mx±{t)) 



±- 



2' 



so 



Hal{t)=} 

where C{t) — ia/2 and 



C{t) \ 
'Cit) j 



(35) 



(36) 



. ^ »fl[A(0-zr/2]-r!fl(t)(A-»r/2) 

" [A(t) - zr/2]2 + r!|(t) ■ ^"^'^ 

Then, the Hamiltonian Ha(t) = Hao + Hai takes the form 

h 
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Hait) 



~A{t) nn{t) + 2C{t) , 
nnit) - 2C(t) A{t) -iT I- 



The practical realization of this Hamiltonian is not 
straightforward. In particular the off-diagonal terms are 
not the complex conjugate of each other unless the real 
part of C(t) becomes zero. We shall explore in the follow- 
ing subsection the possibility to manipulate this result by 
playing with different phases as in Eq. (P7| . 



B. H^{t) applied to a decaying two-level atom 

For the decaying two-level atom, using Eq. (I?f)) with 
phases ^+ = £,+ {t) and ^_ = £,^{t) associated with 
\x+it)) and \x-it)), wc find 



h 

-sin^(f)e+-cos2(f)t 



-cos2(f)e+-sin2(f)i 



(39) 



The phases in the matrix elements H^a.i2{t) and 
H^a.2i{t) only affect the first terms, which are equal. In 
general the manipulation of the phases is not enough to 
make the non-diagonal terms complex conjugate of each 
other since this requires not only Im[(^_ — ^_|-) sin a] = 
but Re[C(i)] = too. We also add potentially complex 
terms in the diagonal that again could complicate the 
physical realization. 

In summary, the phase manipulation does not help to 
implement the shortcut. In some parameter regimes, 
however, an approximation to Hai that leads to essen- 
tially the same results may be easily realized, as discussed 
next. 



C. Forced population inversion 

We study now the forced coherent decay from the up- 
per level of a two-level system with slow spontaneous de- 
cay. This decay may be driven and accelerated adiabati- 
cally with a "rapid" adiabatic passage (RAP) technique. 
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sweeping the laser frequency across resonance. The ad- 
jective "rapid" here could be misleading: it simply means 
"faster than the spontaneous decay" but, as the approach 
is adiabatic, it fails for short enough times. The adia- 
baticity criterion is worked out in the appendix. To go 
beyond the time limits imposed by the breakdown of adi- 
abaticity, shortcut techniques may be applied. 




FIG. 1: Population of the ground state, Pi{t) (solid line), and 
of the excited state, P2{t) (dashed line), for the Hamiltonian 
Haoit). Parameters: P = 27r x 2 MHz, a = (27r)^ x 0.01 GHz^ 
b = (27r)^ X 0.00025 GHz^ and Qo = 2tt x 100 MHz. 

We consider a linearly chirped Gaussian pulse with de- 
tuning A{t) = LjQ — uj{t) = —2bt and Gaussian Rabi fre- 
quency r^H(t) = r2oe~°* • 

The initial conditions are Pi{0) = and P2{0) = 1. In 
Fig. [T]we show that the application of a RAP pulse with 
Hao{t) is only partially successful. Note the slow spon- 
taneous decay before and after the pulse, and a faster 
forced transition during the pulse around t = 0. Since 
the pulse duration is too short, adiabaticity fails. Fig. 
[5] shows the fast full population inversion when adding 
the Hamiltonian Hai{t) in Eq. (155)1 . This Hamiltonian 
has off-diagonal terms with real and imaginary parts de- 
picted in Fig. [21 Whereas the imaginary parts, the bigger 
bumps in Fig. [31 are realizable (2^, the real parts con- 
stitute a non-Hermitian contribution. They are however 
small, and an approximation of Hai{t) neglecting them 
provides essentially the same dynamics, as shown in Fig. 
[51 This remains valid in the strong-driving regime in 
which r -C ilo and the natural lifetime is large compared 
to the duration of the forced decay. 



IV. INVARIANTS BASED INVERSE 
ENGINEERING 

Lewis and Riesenfeld pH proposed the use of dynami- 
cal invariants of a quantum mechanical system to perform 
expansions of arbitrary time-dependent wave functions 
by superposition of eigenstates of the invariant. This 
may be generalized to non-Hermitian Hamiltonians [29l - 
[3lj |. We shall assume that for a Hamiltonian Ho{t) with 
the features described in Sec. II Al there is a generalized 



':o.5 




FIG. 2: (Color online) Population of the ground state, Pi{t) 
(solid line), and of the excited state, P2(t) (dashed line), for 
the total Hamiltonian Ha{t), coinciding with the populations 
Pi(t) (triangles) and P2(t) (circles) when Ha{t) is approxi- 
mated by neglecting Re[C(t)]. Parameters as in Fig. [T] 




FIG. 3: Real (solid line) and imaginary (dashed line) parts of 
C. Parameters as in Fig. [1] 



invariant I{t) that satisfies 



0, 



(40) 



so that ^(*(t)|/(t)|*(i)) = 0. Note that this is not an 
ordinary expectation value {'^ {t)\I {t)\'i {t)) / {"^ {t)\^ [t)) , 
in this sense the concept of generalized invariant differs 
from the one for Hermitian Hamiltonians. 

Let us assume also that I{t) has a non-degenerate 
complete biorthonormal set of instantaneous eigenstates, 
{\4'n{t))^ \i'n{t))}^ whcrc u varics from 1 to N , that sat- 
isfy 

/(Ol^nW) - In{t)\Mt)), (41) 
lHmn{t)) = imMt)), (42) 
{4^m{t)\Mt)) - Smn, (43) 

J2\Mt)){^n{t)\ = 1. (44) 

n 

We can write the general solutions of the Schrodinger 
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equations for HQ{t) and Hl{t), Eqs. ([5]) and as 

|*(t)) = ^d„e-"W|^„(t)), (45) 

n 

\i>{t)) - 5]d„e-"W|^„(t)), (46) 

n 

where the coefficients and {c?„} do not depend on 

time, and the generahzed Lewis-Riesenfeld phases are 



Unit) 



.d_ 

''dt' 



m') 



Mt'))dt'. (47) 



Inverse engineering techniques rely on designing the in- 
variant eigenvectors and phase factors first, possibly tak- 
ing into account partial information on the structure of 
the Hamiltonian, and then deducing the Hamiltonian 
from them. 



V. CLASSICAL PARTICLE IN AN EXPANDING 
HARMONIC TRAP 

It is possible to study a classical particle with position 
q{t) and momentum p(t) in a harmonic trap as a formal 
quantum two-level system with non-Hermitian Hamilto- 
nian, by rewriting the classical canonical equations of 
motion in matrix form 29, 30]. The Hamiltonian of a 
classical harmonic oscillator with a time dependent fre- 
quency w(t) is 



Hho{t) 



2m 



(48) 



where m is the mass of the particle. We shall consider 
an expansion from wq = i^(0) ai t — Q to ujf — uj{tf) at 
the final time t = tf, with Uf < ojq. The corresponding 
classical canonical equations 



= 

Pj = 

can be written as 



dpj m 
dHho 







I/to 
-muj'^it) 



Pit) 



(49) 
(50) 

(51) 



due to their linear dependence on q and p. Multiplying 
both sides of the equality by i we obtain a Schrodinger- 
like equation {H = 1) with the "effective" non-Hermitian 
Hamiltonian 



1/m 



(52) 



This is a useful but formal analogy, since this Hamilto- 
nian does not have units of energy, in fact different matrix 



elements have different units as the "state vector" com- 
ponents q and p have also different units. Nevertheless 
we may apply the generalized invariant theory, and ex- 
pand the state vector in terms of formal eigenvectors of 
the generalized invariants. Defining (soj 



_ / b{t) c{t) 



" 1 -a{t) -b{t) 
and imposing Eq. (j40p . without h, 

UJQ 1 



bit) 
c{t) 



£»2 [t) Wo 



-Q\t) 



-1 

Wo 
^2 



g{t)Q{t), 



WqTO ' 



(53) 

(54) 
(55) 
(56) 



where the dimensionless scaling function git) satisfies the 
auxiliary equation 



Q{t)+^''{i)Q{t) 



(57) 



which is the Ermakov equation, the same equation for the 
scaling function that defines the invariants in the expan- 
sion of the quantum harmonic oscillator For I{t), 
whose eigenvalues are I± = ^i, the eigenstates are 



lV'±(i)) 



cit) 
±i [l±ih{t)] 



(58) 



in the basis used in Eq. ((52|) . The Lewis-Riesenfeld 
phases a±{t) are 



a±{t) - (V^±(i'; 



= I In 




d 



c{t) /•* 1 



dt'. (59) 



Then, the phase-space trajectory is given by 

= d+e-+(*)|V^+(t))+d_e»-(*)|^_(t)) 



= R 



Q{t) 



g{t) cose(t) 
sin 6{t) + rag{t) cos 0{t) 



(60) 



where R ~ 2r c{0) / mujo is a distance, d-^- ^ d*_ = 
rexp(i0o) can be determined by the initial conditions 
at t = 0, and 



d{t) = LUo 



1 



9' it') 



dt' 



(61) 



with ^0 the initial phase. 

Imposing the boundary conditions ^3(0) = 1 and g{0) ~ 
0, and gitf) — (cjo/w/)^/^ and gitf) — 0, which consis- 
tently with the Ermakov equation imply ^i(O) = gitf) = 
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0, we find := E{t = 0) = (jjlR^m/2 and Ef := 
E{tf) = ujfEo/LOQ. In other words, these boundary condi- 
tions guarantee that the value of the classical adiabatic 
invariant E{t)/uj{t) at initial and final times coincides, 
even though it may take different values at intermediate 
times. 

To design the process, Q{t) has to be interpolated at 
intermediate times. We assume here a polynomial form, 
Q{t) ~ Y^n=o ^nt"^ where the coefficients a„ are fixed 
from the boundary conditions. Then we get uj(t) from 
the Ermakov equation. 



vided application examples. Experimental implementa- 
tions are at reach. Related open questions are the ap- 
plication of similar concepts to master equations, or de- 
veloping means to implement arbitrary non-Hermitian 
interactions. Another interesting research avenue is to 
combine shortcut techniques with optimal control [l6l.f20l| 
taking into account physically imposed constraints. 
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g^t) 



m 



(62) 



In Fig. 2] we have represented the shortcut trajectory 
in phase space between the initial and final times, t = 
and t = tf, for the frequency aj(i) given by Eq. (1621) . 
We have also added a period Tq — 2tt/ujo before t = 0, 
and a period Tf — ^Ti/ujf after t = tf, for which the 
particle evolves for fixed and w/, respectively, so as 
to depict complete initial and final ellipses. The shortcut 
trajectory that connects the initial and final ellipses is 
clearly not an adiabatic path, that would be formed by 
a succession of slowly varying ellipses from the initial to 
the final one. 
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Appendix A: Adiabaticity condition for 
time-dependent non-Hermitian Hamiltonians 
applied to a decaying two-level atom 




FIG. 4: (Color online) Parametric velocity-position trajec- 
tory. The initial ellipse (solid blue thin line) and the final 
ellipse (red dashed line) are connected by the shortcut trajec- 
tory (green solid thick line). Parameters: = Q correspond- 
ing to go = q[t — Q) — 1 fim and vo — v{t = 0) = /^m/ms, 
cjo = 27r X 250 Hz, ujf = 2tv x 2.5 Hz, tf = 25 ms, and the 
mass of an atom of Rubidium-87, m — 1.44 x 10~^^ kg. 



The adiabaticity condition for time-dependent Hermi- 
tian Hamiltonians is given by 



\{n{t)\dtm{t))\ < ^\E.n{t) - E„^{t)\, n ^ m, (Al) 



in terms of instantaneous eigenstates and eigenvalues. 
Following closely its derivation in [s^l we generalize it 
for non-Hermitian Hamiltonians as 



\{n{t)\dtm{t))\ < ^\Enit) - Emit)\. (A2) 



For the two-level decaying atom this condition is 



\{n+it)\dtn^m <$: -\E+{t) - E^it)\. 



(A3) 



Introducing here Eqs. (f^ . (15(11) and ([55t . the adiabatic- 
ity condition for this system takes the form 



VI. DISCUSSION AND CONCLUSION 



2\na{t)\ < \n{t)\, 



(A4) 



We have generalized shortcut to adiabaticity tech- 
niques for non-Hermitian Hamiltonian systems and pro- 



where n{t) = ^ -[T + 2iA{t)]^ + Aflj^it) and flait) = 
-d/2. 
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